Unplugging the Universe: the neglected electromagnetic consequence of decoupling 
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This letter concentrates on the non-equilibrium evolution of magnetic field structures at the onset 
of recombination, when the charged particle current densities decay as neutrals are formed. 

We consider the effect that a decaying magnetic flux has on the acceleration of particles via the 
transient induced electric field. Since the residual charged-particle number density is small as a 
result of decoupling, we shall consider the magnetic and electric fields essentially to be imposed, 
neglecting the feedback from any minority accelerated population. 

We find that the electromagnetic treatment of this phase transition can produce energetic electrons 
scattered throughout the Universe. Such particles could have a significant effect on cosmic evolution 
in several ways: (i) their presence could delay the effective end of the recombination era; (ii) they 
could give rise to plasma concentrations that could enhance early gravitational collapse of matter by 
opposing cosmic expansion to a greater degree than neutral matter could; (iii) they could continue 
to be accelerated, and become the seed for reionisation at the later epoch z « 10. 

PACS numbers: 98.80.-k, 98.80.Jk, 91.25.Cw, 94.20.wc, 98. 58. Ay, 98.62.En, 98.70.Sa 



In this article, we consider the effect that a decaying 
magnetic flux has on the acceleration of particles via the 
transient induced electric field. We do not address the 
origin of the magnetic field, but merely assume that it 
exists (see [l[ and the references therein). 

The Universe underwent recombination at redshift z ~ 
1100, in which the majority of charged particles (i.e. the 
current-carriers) disappeared as neutrals were formed. 
This can be understood in the context of simple cosmo- 
logical models (e.g. @, H, 01). In general, such mod- 
els do not attempt to describe in detail the actual non- 
equilibrium processes themselves, but instead simply ac- 
commodate their broad consequences. 

However there is always a risk, when presenting a series 
of (quasi) equilibrium states as an evolutionary sequence, 
of overlooking some essential aspect of the physics that 
permits the very transition between states. In this arti- 
cle we develop the consequences of describing the tran- 
sient physics that occurs at decoupling: current densi- 
ties that decay must inevitably cause an induced elec- 
tric field, which must accelerate any remaining ambient 
charged particles. We shall show that this physics can- 
not be contained within the fluid plasma models that 
are conventionally used to describe the transition; this 
is why such essential physics is missing from the overall 
understanding of universal evolution. 

Electromagnetism versus fluid models. — Electromag- 
netism tells us that magnetic fields are sustained by dis- 
placement or charged particle currents. Prior to recom- 
bination the high conductivity of the plasma makes the 
latter the most plausible source of magnetic curvature 
[B|, After decoupling only a small fraction of elec- 
trons x e ~ 10~ 5 [7| persists, meaning that the conti- 
nuity of the total current must inevitably involve the 
displacement current replacing the particle current, and 
consequently the evolution of the pre-decoupling mag- 
netic structures. The accompanying induced electric field 



will accelerate any remnant charged particles, providing 
a minority source of energetic plasma that may well be 
significant in the onward evolution of the cosmos. 

Suppose that a typical magnetic field Bq pervades a re- 
gion of characteristic dimension Lq and that the particle 
current which is sustaining it vanishes. How does the sys- 
tem evolve from this point onwards? Electromagnetism 
tells us that if the particle currents are interrupted, then 
the overall system will react with an induced electric field 
in an attempt to oppose the change. The size of the in- 
duced electric field is related to the timescale over which 
the magnetic flux is changing. The subtlety here is en- 
suring that cosmological models permit the appropriate 
timescales in the evolutionary processes. Fully electro- 
magnetic descriptions are entirely appropriate, but those 
that depend on magnetohydrodynamical (MHD) plas- 
mas may well be unable to accommodate the necessary 
physics at transition, and as a result, underestimate the 
significance of the decoupling itself. 

In the fully electromagnetic description, the electric 
and magnetic fields E, B and plasma current density J 
are governed by 

VxE = -B, V x B = [i a J + E/c 2 , (1) 

where ' denotes d/dt. Assuming that the current density 
J is negligible as a result of decoupling leads to V x B 
being determined by E/c 2 . The dimensional analysis of 
yields a timescale for changes in the field evolution of 

n = Lq/c (2) 

Hence the magnetic field reconfigures by changes which 
propagate at light speed. 

However, if the plasma is modelled using MHD such 
rapid communication is not permitted, since the displace- 
ment current is omitted from MHD dynamics Hence 
instead of the full electromagnetic equations, the fluid 



plasma of resistivity r\ is governed by 

V xE= B, V xB = fi J, (3) 
E + u x B = r\ J, (4) 

in which u is the plasma velocity. Note that now we can- 
not ignore J in favour of the displacement current, since 
the former is the only way in which V x B is sustained. 
It is clear at this stage that the sudden loss of electri- 
cal current as neutrals form is beyond the compass of 
such a framework; the magnetic field can only communi- 
cate changes at the resistive or fluid-dynamic timescales. 
This can be seen from the form of Ohm's Law in Eq. ([4]) . 
Combining Eq.s ([3]|4]) we obtain 

B = V x (it x B) - fjV 2 B (5) 

in which fj = T)/[1,q = constant. In magnitude terms 

B q /t : u B /L : f)B /Ll or 1 : u t/L q : fjr/Ll (6) 

which leads to a typical resistive time scale Ti given by 

T2 = Lim (7) 

assumed to be the characteristic of resistive MHD; the 
dynamo timescale 

r 3 = L Q /u (8) 

is the only possibility for an ideal (perfectly conducting) 
MHD plasma. 

The induced electric field. — The magnitudes of the in- 
duced electric fields for the different scenarios are given 
by the induction law, which relates the induced electro- 
motive force to the rate of change of magnetic flux $: 

f E - dl = -t = l-tJ B - dA < 9 > 

in which the closed-loop integral is around the path en- 
closing the area A, in the usual way. We have used par- 
tial derivatives in time on the right-hand side because we 
will consider the field at a point in space, rather than 
the motional emf. If we assume a circular area of radius 
Lq through which the magnetic flux is linked, then we 
can identify the typical size of induced electric field by 
dimensional analysis: 

E ~ BL /t d (10) 

where td characterises the time over which the magnetic 
field decays at that location, that is, the timescale for 
the particle current density J to disappear. If we use the 
timescales associated with the electromagnetic and MHD 
models, then we arrive at the following estimates of the 
induced electric field: 

electromagnetism E\ — cBq (11) 
resistive MHD E 2 = L /t 2 B = f)/L a B a (12) 
ideal MHD E 3 = u B (13) 



Comparing the magnitudes of the fields by forming E\ : 
E% : E3 we have 

c:B : fj/L B : u B or 1 :f]/L a c :u /c (14) 

Since MHD (ideal or resistive) must have uq <C c, it is 
clear from this that the induced electric field permitted 
by the fully electromagnetic model is the most significant 

The key issue here is that the size of the induced elec- 
tric field is determined by the rate of change of the mag- 
netic flux. Only the electromagnetic model can accom- 
modate the rapidly varying current density that results 
from recombination (a quantum, rather than fluid, ef- 
fect). The fluid plasma models tie the evolution of the 
magnetic field explicitly and directly to the fluid motion 
via u. This is particularly acute in ideal MHD, where 
the flux linked by a fluid element can never change, and 
so the only physical process that allows flux to decay 
overall is cosmic expansion. There is a subtle point here 
in ideal MHD: the frozen- in flux condition means that 
the magnetic flux at a point can change only at the ex- 
pense of plasma motion; however, a decrease at one point 
leads to an increase somewhere else as the fluid rarifics 
and compresses. Hence ideal MHD cannot in fact ac- 
commodate recombination, since the latter causes an un- 
avoidable overall loss of magnetic flux when the current 
density (essentially) vanishes - it is not simply a case of 
rearrangement . 

Moreover, the fluid model is further compromised by 
the fact that the u that appears is the plasma fluid veloc- 
ity: the very species that is being eliminated as a result of 
the recombination. Hence tying the induction field to a 
vanishing fluid species cannot recover the correct physics 
at decoupling. This explains why the electric field pro- 
duced in such models is deemed wholly negligible. 

The consequence of electric induction. — We will now 
explore the consequences of an induced electric field, con- 
sidering only the electromagnetic case. It is vital to dis- 
tinguish between timescales, however, in order to avoid 
confusion. 

The timescale for transmission of field changes (in- 
cluding boundary conditions) is the electromagnetic one, 
given by t\ in Eq. {2]); we will drop the subscript 1 from 
now on. 

There is also the characteristic time for the decay of 
the current density, m, which could be very short, since 
it is governed by atomic processes: for example, by the 
rate coefficient for the attachment of free electrons to 
positive ions. The process is strongly density dependent, 
and therefore influenced by local conditions. 

Finally, we can identify a timescale tr over which re- 
combination is completed across the universe. This is 
taken to be of the order of 100 kyr ~ 10 12 s. 0] 

Hence at a location, changes in the current density due 
to neutral formation (evolving over the timescale td) are 
communicated throughout that location by electromag- 
netic signals, which take a characteristic time r to reach 
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all parts, tr is the timescale for all such locations in the 
universe to complete this process. 

A simple example: particle acceleration. — Dimen- 
sional analysis of Eq. ((9]) gives the electric field magni- 
tude in terms of the magnetic field and the characteristic 
length and time: 

E = BL /t d (15) 

In order to quantify the size of the induced electric field, 
we must supply typical values for B, Lq and r n ■ A typ- 
ical magnetic field value is B ~ lCr 12 T 0, \U^. The 
scale-length is more challenging: this could range from 
10 11 m (an astronomical unit) to 10 16 m (a parsec), with 
the latter being more realistic in the context of structure 
boundaries visible in [jj [l(| • 

The characteristic time for the recombination to take 
place at a point is not the recombination width tr, but 
the timescale that arises self-consistently from the rate 
equation for the formation of neutrals. In general terms, 
the rate equation governing the conversion of one species 
into another via collisions with electrons takes the form 
h a = anbn e where here n a plays the role of neutral num- 
ber density, nj, the positive ion number density, and n e 
the electron number density. The rate coefficient a for 
recombination 11, 12. [H is typically 10~ 16 m 3 s -1 , yield- 
ing 

t d ~ (na)- 1 (16) 

where n is a typical number density at the recombination 
epoch. Taking n ~ 10 8 m~ 3 yields td ~ 10 8 s. 

Suppose we have a uniform axial magnetic field B(x) 
within some radius Lq which is supported by an edge- 
located current structure which vanishes at the time of 
recombination. (A helpful analogy would be a magnetic 
solenoid.) The magnetic field is given by B(x) = zB 
and VxB = except where the current is carried. When 
the current density is switched off, the induced electric 
field is non-zero everywhere (except perhaps within the 
original current carrying region) as can be seen from 
Eq. ([1]). We shall assume that the induced field prop- 
agates into the uniform magnetic field region at speed c, 
and our analysis will be concerned with the evolution of 
the system once the current collapse has been communi- 
cated across Lq. This is a reasonable assumption, since 
we wish to follow the energetic evolution of an inertial 
test particle initially at rest when the current collapses. 
Of course, any candidate for acceleration by such fields 
will initially be in equilibrium with the CMBR, but this 
straightforward case will highlight the essential physics. 

Let us assume that B retains a predominantly ax- 
ial character during the acceleration of the test particle; 
there will of course be evolving magnetic curvature, but 
in the spirit of keeping the problem simple, we will set 
this aside. Hence the induced electric field will remain 
primarily azimuthal. Therefore 

p - l {d/d p ){ P E^)«-B {t) (17) 



where (p, <j>, z) are standard cylindrical coordinates. The 
right-hand side of Eq. (|17p is independent of p and so we 
can integrate immediately to get 



Et(p, t)w-So(«) P/2. 



(18) 



Given these approximate fields, we can now turn to the 
trajectory of an electron that is subject to them. It is ap- 
propriate to consider the single-particle behaviour here, 
since we are describing the evolution of a minority pop- 
ulation of charged particles that escape elimination from 
neutral recombination. We know that the ionization frac- 
tion drops by 5 orders of magnitude at recombination 
0; it is reasonable then to consider single-particle ef- 
fects without taking into account any bulk feedback from 
the motion of such charges. Note that the particle being 
accelerated here needn't necessarily be an original mem- 
ber of the particle current ensemble that sustained the 
pre-decoupling magnetic field. 

The relativistic Lorentz equation for an electron mov- 
ing under imposed electric and magnetic fields is 



dp/dt = -e(E + v x B) 



(19) 



where p = jm e v is the electron 3-momentum, and 
7 = I/a/I — v 2 /c 2 . The reference (or lab) frame here 
is the one in which the magnetic field is predominantly 
axial, and the electric field azimuthal, consistent with our 
model configuration. 

If this particle acquires the drift speed Ud = 
E x B / B 2 in the lab frame then the Lorentz force van- 
ishes, and the particle is free, having energy £ in addi- 
tion to its rest energy given by £ — (7 — l)m e c 2 . We 
can take £ to be the maximum energy gain for a parti- 
cle initially at rest (or in thermal equilibrium at decou- 
pling) before being subjected to the induced electric field 
that must accompany the collapsing current density. For 
E/B «5x 10 7 ms _1 , then 7 w 1.014 and so such an elec- 
tron has gained approximately 1.4 % of its rest energy, or 
about 7 keV, as a direct result of the decoupling process. 

This result can be achieved in a different way. Consider 
a thermal electron at decoupling, moving in a Larmor 
orbit in our typical magnetic field of 1 pT. Using vj_ « 
3 x 10 5 ms" 1 and to c = eB§jm e « 0.1 Hz, the Larmor ra- 
dius Rl — v/u c at decoupling is Rl ~3x 10 6 m. Since 
the induced electric field is predominantly azimuthal, 
the change in azimuthal speed 6v resulting from the az- 
imuthal acceleration is Sv ~ (e/m e )ESt, leading to an 
increase in Larmor radius of 5Lr ss (e/u) c m e )E8t. Tak- 
ing the outward drift speed of the particle as the rate of 
change of its Larmor radius, we arrive at di?i/dt ks E/B, 
as before. Substituting for E using Eq. (fT8|) , we can solve 
for the time evolution of the outward drift: 



R L {t)~Rw [B /B{t)) 1/2 



(20) 



where -Bo is the magnetic field just prior to decoupling. 
If we assume that the current density decays exponen- 
tially with characteristic time td, then it is clear that 
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the Larmor radius of the particle reaches lpc after a time 
t ~ 40t£), at which point it is travelling with a radial 
speed di?i/di w Rl /(2td) ~ c/6, as before. Note that 
this outward drift speed dominates the azimuthal com- 
ponent. 

This modelling underlines the fact that the elimina- 
tion of magnetic energy at decoupling creates an elec- 
tromagnetic impulse which accelerates ambient charged 
particles. 

Conclusions. — In this letter we model the transient 
magnetic and electric fields at the time of recombination, 
distinguishing between three critical timescales: (i) r, the 
time taken for fields to accommodate the evolving current 
density; (ii) td, the time for recombination to take place 
in a localised region; and (iii) tr the time window for 
decoupling to be completed across the universe. 

Dimensional analysis shows that only a full electromag- 
netic model allows r to be consistent with the appropriate 
physics: recombination is a relatively rapid event driven 
by the relevant rate coefficients and number densities of 
participating species, and not by a global fluid (MHD) 
model that cannot incorporate rapid changes in current 
density, and in the case of ideal MHD, cannot allow a 
global change in the magnetic flux. Such MHD mod- 
els inevitably link field evolution to bulk fluid motion, 
thereby imposing an unacceptabe evolution timescale. 

The critical features of our model can be summarised 
as follows. A minority cohort of charged particles (here 
taken as electrons), initially in thermal equilibrium with 
the matter in the Universe, survive the recombination 
process at the boundary between magnetic domains. The 
rapidly changing current density caused by the elimina- 
tion of the majority population of free charged particles 
(recombination) induces an electric field that accelerates 
any that are left over. The maximum energy that such 
particles can acquire in this process occurs when they 
move at the drift velocity, corresponding to an energy of 
approximately 1.4% of their rest mass. 

We have identified a process in which a population of 
energetic particles is produced as a direct result of de- 
coupling. Such particles, although very much a minority 
fraction of the neutral number density, are clearly not in 
thermal equilibrium with that majority; this may leave 
a signature in the CMBR, albeit at very low frequency. 
Moreover they are sufficiently energetic to cause an ele- 
ment of reionisation of neutral matter in the decoupling 
era, should the collisional environment be appropriate, 
and in so doing, perhaps delaying the end of the recom- 
bination era. 

More significantly we have described a mechanism 
which can seed the Universe with high-energy charged 
particles from a single, localised decoupling event. Such 
particles could in principle be accelerated to even greater 
energies by further encounters with induced fields at dif- 
ferent locations, since it is clear that tr 3> td. Hence 



our model may offer a route whereby cosmic ray particles 
could be created. 

Finally, although in a minority, the persistence of such 
energetic plasma beyond the decoupling era may acceler- 
ate the condensation of matter into early structure forma- 
tion, since long-range Coulomb forces may be sufficient 
to bias gravitational self-collapse. 

A full numerical simulation of the magnetic field evolu- 
tion, energetic particle production and the observational 
consequences is currently underway [14j and will be re- 
ported later. 
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